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Abstract. In this note, we discuss the quantum Hele-Shaw flow, a random 
measure process in the complex plane introduced by the physicists P.Wieg- 
mann, A. Zabrodin, et al. This process arises in the theory of electronic 
droplets confined to a plane under a strong magnetic field, as well as in the 
theory of random normal matrices. We extend a result of Elbau and Felder 
|B] to general external field potentials, and also show that if the potential is 
C 2 -smooth, then the quantum Hele-Shaw flow converges, under appropriate 
scaling, to the classical (weighted) Hele-Shaw flow, which can be modeled in 
terms of an obstacle problem. 



1. Introduction 

A weighted distribution of TV-tuples of complex numbers. The intention of this 
paper is to provide a rigorous mathematical treatment of the theoretical physics 
model due to Wiegmann, Zabrodin, et al., which is concerned with the behavior 
of a finite number of electrons (charged fermions) in a strong magnetic field in the 
context of Quantum Theory. In part, this was done very recently by Elbau and 
Felder. We will say more about this below. 

Let N be a positive integer, and introduce the notation 

z [i,N] = (zi,z 2 ,z 3 ,...,z N ) 

for a point in iV-dimensional complex Euclidean space C^. As a real vector space, 
we have the identification C N ~ R 2N . We let X 2 n be the standard Lebesgue 
measure on C^, so that in particular, 

dA 2 jv(z[i,jv]) = dA 2 (zi) • • ■d\ 2 (z N ). 

The van der Monde determinant is the quantity 

A(z [1 , A r ] ) = Y[ ( Z 3 ~ Z >=)> 
j,k: j<k 

where we tacitly assume that the parameters j and k range over the set {1, . . . , A''}. 
Let $ : C — > M be a C 2 -smooth function, which grows so quickly that 

(1.1) $(z) > A \og\z\ +0(1) 
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holds for all positive real constants A (with O(l) dependent on A). Let j3 denote 
a positive real parameter. We let /xg^ be the finite positive Borel measure on C 
given by 

(1.2) dfip^iz) = exp{ -/3$(z)}dA 2 (z), z e C, 
We shall consider the probability measure 

1 2 

(1.3) dnAT,/3,$(zn jv]) = — A(z 1)iV ]) d/i /3 .$(zi) • • -<lu,p^{z N ), 
where Zn is the normalizing constant: 

(1.4) Z N = |A(z ljAr i)| 2 d/z / 3 i $(zi) • • •d/i ( 3,*(2Ar). 

We intend to study the typical behavior of (|1.3|) for large values of N. To make 
this more precise, let n be a positive integer with 1 < n < N, and consider the 
marginal distribution of 

Z[l,n] = (z 1 ,Z 2 ,Z 3 ,...,Z n ), 

as N — ► +oo while n is kept fixed. This n-point marginal distribution is given by 
the probability measure 

(1.5) dn^g )# (z [1)n ]) = -J- I ^ |A(z [liA r])| 2 d/i /3 ,$(z„ + i) ■ • ■ djU/s,* (zjv)j 

x dfj,p t &(zi) ■ ■•d/i l g,$(z n ). 

The measure Tfop $ obtained by multiplying Il^g $ by Nl/(N — n)\ is called 
the n-point correlation measure. The distribution l|1.3f) appears from at least two 
different contexts: one is the theory of eigenvalues of random normal matrices, and 
the other a physical configuration of N electrons localized to a two-dimensional 
plane, while being exposed to a strong magnetic field (the strength is regulated by 
the parameter 0) perpendicular to the given plane. In the second instance, /3$ is a 
scalar magnetic potential, z\,...,Zn are the locations of the N electrons, and one 
looks for time-independent solutions to the Schrodinger equation. The resulting 
cloud of electrons is an instance of the Aharonov-Bohm effect, where the magnetic 
potential rather than the magnetic field manifests itself. For details, we ask the 
reader to consult the work of Wiegmann, Zabrodin, et al. P, El' CB1> an d 
|21| . In the first instance (random normal matrices), we should mention first that 
normal N x N matrices have the decomposition M — U*DU, where D is diagonal 
(consisting of the eigenvalues), and U is unitary, that is, U*U = UU* = I, where 
I is the identity matrix. To focus on the eigenvalues, then, one should remove as 
much as possible of the unitary part; in other words, we should mod out with 
respect to U [N)Jli d(N), where U (N) is the group of unitary N x N matrices, 
and U a{N) is the subgroup of unitary diagonal matrices. If this is done carefully, 
and a weight is introduced to get a probability measure on the normal matrices, 
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we get a joint distribution of the eigenvalues of the type l|1.3|) . For details of this, 
see [B] and @|; see also Mehta's book |T5) . 

The standard models for the eigenvalue distribution of self-adjoint matrices 
are by now quite well-known (see |18|. |13p. In the limit, we get such models from 
our probability measure (|1.3|l by letting <&(z) tend to +oo on C \ 1. For instance, 
the semicircle law of Wigner can be obtained by taking such a limit with the 
appropriate interpretation. 

Quantum Hele-Shaw flow vs Hele-Shaw flow. We show that the n-point correlation 
measure 

(n) Nl n (n) 

is given in terms of a determinantal process involving reproducing kernels of poly- 
nomial subspaces of a Bargmann-Fock type space, and that it defines a monotonic 
growth process in N (for fixed /3 and $): 

p(") < p(«) ivr _ i 9 o 

L N,{3,<S> — 1 JV ~~ 1 i z > °> 

For this reason, and for reasons which will be made clearer below, we call this 
Quantum Hele-Shaw Row. The analogy with classical Hele-Shaw flow was made 
earlier by Wiegmann and Zabrodin. The ratio N / (3 seems like a natural growth 
parameter for the process. For mathematical reasons, however, it is more natural 
to use 

, s N - 1 

(1-6) r= — 

instead. 

The weighted energy of a compactly supported Borel probability measure a 
is (for a given positive parameter r) 

£*M (7) =5; ^{2rlog^^+*(*)+*(iD)|d£r(2)d<7(ti;); 

the unique compactly supported Borel probability measure which minimizes this 
energy is called the extremal measure, or equilibrium measure, and denoted by 
er$. T , or simply by er, if no misunderstanding is possible. The existence and unique- 
ness of the extremal measure are treated in Chapter 1 of Saff's and Totik's book 
d 

Our first main result reads as follows. 

Theorem 1.1. If t, < t < +00, is kept fixed, and (3, N are connected via |J.6|) . 
then, as N — > +00, the marginal distribution has the weak- star limit 

dn jV,/3,<I.( Z [l,n]) &r( z l) ■ ■ ■ <&( z n)- 



This theorem is essentially as in Johansson's paper |13j . but with important 
additional technical details. Elbau and Felder |S] were the first to realize that 
Johansson's argument extends to the case of normal random matrices. The present 
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version is more general - Elbau and Felder required that the potential be infinite 
in a neighborhood of infinity. 

We discuss the extremal measure in Section |3 and connect it to an obstacle 
problem in Section^ This leads to our second main theorem. For the formulation, 
we need the following standard notation: A is the standard Laplacian in the plane, 
while given a subset E of the complex plane C, 1e is the characteristic function 
of the set E. 

Theorem 1.2. For each t, < r < +oo, there exists an unbounded open set 2?$(r) 
in C whose complement is compact, such that the extremal measure has the form 

da^ T (z) = — l C \x.*(r)(z) A$(z) dA 2 (z), z e C. 

Moreover, as r increases, the open sets 2?$ (r) form a continuous decreasing chain, 
and the measures tS$ iT grow with r. 

This theorem justifies the claim of Wiegmann, Zabrodin, et al., that in the 
limit one gets a monotone family of domains (and not only measures) in the 
complex plane. The proof is an application of the general theory of free boundary 
problems due to Caffarelli et al. 

We should add that, if the weight function $ is real analytic, then the bound- 
ary 92?$ (r) consists of real analytic curves, except for finitely many cusps (pointing 
outward) or so-called contact points. The evolution of 2?$(r) is quite similar to 
the Hele-Shaw flow on hyperbolic or weakly hyperbolic surfaces described in [T2*] 
and (TT]. In fact, the present situation corresponds to Hele-Shaw flow on surfaces 
where the metric may be negative in portions of space, since $ is a kind of metric 
potential, with associated metric 

ds(z) 2 = A$(z) |dz| 2 . 

Generally speaking, Hele-Shaw flow models the expansion of viscous fluid in a 
surrounding medium (which we may think of as vacuum), as new fluid is injected 
at a constant rate at one or several source points. In the present setting, we have 
a single source point at infinity. 

A final remark we wish to make here is that essentially all results of the 
paper remain valid with little change if the C 2 -smoothness of $ is weakened to 
local H /2,p -smoothness, for some p, 1 < p < +oo. 

2. Quantum Hele-Shaw flow and Bargmann-Fock type spaces 

Weighted Bargmann-Fock spaces and reproducing kernels. Let /i be a positive 
finite Borel measure on the complex plane C, with 

0< / |z| 2n d/i(z) < +oo, n = 0,l,2,3,.... 
Jc 

The Hilbert space L 2 (C,/i) consists of the square-integrable complex- valued func- 
tions with respect to /x, where two functions are identified if the coincide /j-almost 
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everywhere. We shall assume that fj, is supported on an infinite set, so that £ 2 (C, u) 
becomes infinite-dimensional. The above integrability condition ensures that all 
(analytic) polynomials belong to L 2 (C, u). Let P 2 (C, n) denote the closure of the 
linear space of polynomials in L 2 (C, /z). Also, for N = 1,2,3,..., let P^(C, fi) be 
the iV-dimensional linear space of polynomials of degree < N — 1, supplied with 
the Hilbert space structure of L 2 (C, n). The point evaluations in C are linear func- 
tionals on P^(C,fi), and thus given by the inner product with an element of the 
space: 

p(z) = (p,K^(-, z)) L 2 { ^ = K fljN (w,z)p{w)dfj,(w), zeC, 

Jc 

where, for fixed z € C, the function K^^i^z) is an element of P^(C, /i), and so 
is p. The function K^.n is called the reproducing kernel for the space P^(C,/z). 
The monomials 

(2.1) e j (z)=z j -\ j = l,2,...,N, 

form a vector space basis in P^(C,yLt), but they need not be orthogonal to one 
another. To rectify this, we may apply the Gram-Schmidt process: let ipi = e%, 
and for k = 2, 3, . . . , N, put 

fc-i 

(2.2) cpk = ek + Aj.fc ej, 

where the complex scalars Aj,fc are chosen so that (fk becomes perpendicular to all 
the vectors ex, . . . , e^-i- The functions ipi, if2, y?3, ■ ■ ■ obtained in this fashion 
form an orthogonal sequence. To get an orthonormal sequence, we normalize: 

^)= /l| (Z) ' zeC - 

\m\\L 2 (ji) 

The reproducing kernel may then be conveniently expressed in terms of these 
orthogonal functions: 

N 

(2.3) K^ N (z,w) = ^2c/> j (z)j) j (w), (z,w) eCxC. 

Matrix-valued reproducing kernels and correlation measures. For two vectors 

z [l,n] = (Zl,--->«n) and W[ 1;n ] = (u>i,..., w n ), 
we form the n x n matrix 

( 2 - 4 ) *V,Ar( Z [l,n]: W [l,n]) = [ K H,N (zj , W k )] " fc=1 , 

which we may think of as a matrix-valued reproducing kernel. 

We return to the joint distribution function l|1.3|) . in the slightly altered form 

1 ,2 

(2.5) dll^jv^z^^]) = — |A(z [liJV ])| dfj,{zi)---du.(z N ), 
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where zn,jv] = ■ ■ ■ > z n), and Zn is the normalizing constant: 
(2.6) Z N = |A(z[ liJV ])| 2 d/i(zi) • • -du,(z N ). 



We fix an integer n, 1 < n < N, and split Z[i j\n = ( z [i,n], z ]n,N]), with Z[i, n ] — 
(zi, . . . , z n ) and Z] n ,N] = ( z n+i, ■ ■ ■ j z n)- The n-point marginal distribution mea- 
sure is 

(2.7) dn^(z [ln] ) = — \ j \A(z [hN] )\ 2 d^(z n+1 )- ■ -dn{z N )\ 

x dn(z\) ■ ■ -du.(z n ), 

while the n-point correlation measure is 

(2.8) dr^( Z[1 , n] ) = dn^(z [1>n] ), Z[1 , n] e c\ 

where the integration takes place with respect to the variables z n+ \, . . . , Zn only. 
For n = 1, the measure describes the joint density of the eigenvalues. For 

measures \i of the form l|1.2fl . with A<I> a positive constant in a large disk cen- 
tered at the origin, T^ 1 ^ is rather similar to the majorization function of Aleman, 
Richter, and Sundberg [2]; one major difference, though, is that here, we consider 
backward shift invariant subspaces (P N (C, //)), while in [5], the object of study are 
the forward shift invariant subspaces. 



Proposition 2.1. We have that 
det 



K M .iv( z [i,iV],Z[i,jv]) = -z- \A{Z[1.N]) 



Proof. We first apply simple row or column operations on the representation of the 
van der Monde determinant as the determinant of the N x N matrix with entries 
ej(zk), for j, k = 1, . . . , N. The result is that the van der Monde determinant may 
be calculated based on the entries ipj(zk) instead, whence the result follows from 
the formula defining the determinant in terms of permutations, as in |18| . pp. 
89-95. □ 

The following "linear algebra" description of T^ N in terms of the determinant 
of a matrix reproducing kernel is essentially known (see |18| for the Gaussian 
situation) . 

Theorem 2.2. We have that 

d/i(zi) • ■•d/i(z n ). 



dT %( z [i.n]) = det K^ N (z [1>n] ,z [1>n] ) 
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Proof. We are to show that 

A ( z [i,A'])| 2 dAi(2n+i) ■ -'dfi(z N ) = Z N det [K™ jv( z [i,n]! z [i,n])] ■ 

We follow the approach of |T2], pp. 89-95, and use that by Proposition l2.il the 
identity is valid for n = N , which we take as our starting point. Then succes- 
sive applications of Theorem 5.2.1 in |18j (suitable modified by including complex 
conjugates where needed) yields the result for general n = 1, . . . , N. The proof is 
complete. □ 

(n) 

Corollary 2.3. The n-point correlation measure dT^ N increases monotonically 
with increasing N. Its total mass equals N\/(N ~ n)\. 

Proof. That the total mass of drj^ is N\/{N — n)\ follows from the definition 
(|2.6(l of the normalizing constant Zn- In view of 1)2. 3)1 . we have 

K^n+x{z, w) = K^n(z, w) + (f> N+ i(z) 4>n+i(w), 

so that 

K M,Af+l( Z [l,«] J w [l,n]) = K ^Ar( z [l,n], ~W[l,n]) + *p,JV+l i z [l,n] , ™[l,n] ) J 

where 

* A1 ,jv+i(z [ i ! „],w [l! „ ] ) = ((^N+xizj) (f> N+1 (w k ))" k=1 . 

The matrix $ Ali 7v+i(z[i.n] I z [i.ri]) is positive definite, and so is K™ N (z[i n ^,Z[i n ^). 
As the determinant operation is monotonic with respect to increasing positive 
definiteness (see Proposition 12.41 below), the assertion follows from Theorem 12.21 
The proof is complete. □ 

A linear algebra result. The following proposition is known; for instance, it should 
be possible to derive it from the Horn inequalities, which describe the possible 
eigenvalues of the sum of two Hermitian matrices in terms of the eigenvalues of 
the summands (see, e. g., the survey paper |7j). However, we find it convenient to 
supply instead a short direct proof. 

Proposition 2.4. Let A and B be two positive semi-definite (Hermitian) N x N 
matrices. Then 

max { det (A), det (B)} <det{A + B). 



(N - n)\ J c 



Proof. Clearly, by symmetry, it suffices to show that 

det (A) < det(A + B). 

If det (A) = 0, the assertion follows trivially, as the determinant of a positive semi- 
definite matrix is always > 0. In the non-singular case det (A) > 0, we note that 
A~ x l 2 is a well-defined positive semi-definite matrix, and that A" 1 / 2 B A~ x / 2 is 
positive semi-definite as well. We write 

det (A + B)= det(A) det (/ + A' 1 / 2 B A- 1 / 2 ) , 
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where / denotes the N x N identity matrix. The assertion now follows from the 
fact that the eigenvalues of A~ x / 2 B A^ 1 / 2 are all positive, so that the eigenvalues 
of I + A~ x l' 2 B A~ x l 2 are all greater than 1 (after all, the determinant equals the 
product of the eigenvalues). □ 



3. Weighted Fekete points and the extremal measure 

Weighted Fekete points. Let us analyze the probability measure (|1.3|) . written in 
the form 

1 f N 1 2 

(3.1) dL T Ar i( 3 i $(z[ li j V ]) = 7^— exp \ - $(zj) Ua(z [1iJV ])| dA2jv( z [i,JV]), 

N ^ 3=1 ' 

where z» m = [z\, . . . , zjv), and is the normalization constant. The most likely 
configuration zn m is the one that maximizes 

N *. 

,A], 



exp{-/3^$fe)]|A( Z[li 



among all vectors in C . Considering that 



A' 



Y, + *(**)] = (N - 1) Y 

3,k:j<k 3=1 

we may rewrite this expression: 

(3.2) expj -^^J |A(z M] )| 2 

= n (^(-^-[^o+^^k-^p 

j,k:j<k { V 
For positive real 6*, let E<z>(z, w; 9) denote the function 

(3.3) E 9 (z,w;9) = \z-w\ exp ( - 6 [$(z) + 
we are interested in maximizing 

1 ; E^iz^z^e), 



j,k:j<k 

with 



/3 



2(N - 1)' 

Let M<s> } n{0) be this maximum raised to the power 2/(N(N — 1)): 

, \ 2/(AT(JV-l)) 

Af* lJV (0) = sup J TT E*(zj,z k ;0) \ 

-U.*l6 CW li,fcTj<* J 
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The points Z[i,jv] = z *i n] which achieve this maximum are called weighted Fekete 
points (see |19p. In order that the maximum be assumed, we should add an as- 
sumption on the weight <f>: 

(3.4) E(z,w;9)^0 as max{|,z|, \w\} +00. 

An easy argument shows that (|3.4|) follows from Ijl.lf) with an appropriate choice 
of the parameter A. It follows from Theorem 1.1 p. 143] that M$,^(8) is 
decreasing in N, provided that <& is fixed and 9 is held constant. Moreover, under 
the same conditions, M^, i ^(9) — ► M$(0) as N — ► +00, where 

(3.5) M*(0) =exp(-K$(6»)), 

and k$((9) is results from a certain minimization problem, outlined below. The 
quantity M$(9) can be thought of as a weighted capacity (of the whole plane). 

The extremal measure. Let V c (£) denote the convex "body" of all compactly 
supported Borel probability measure a on C. For a € "P C (C), we introduce the 
energy functional 

(3.6) £Jo-](0) =20 / $(z)dcr(z)+ // log , 1 dcr(z)da(w) 

Jc JJc 2 \z~w\ 

log , 1 , + e $(z) + e $H 1 do-(z)dcr(u)) 

C 2 I \Z-W\ J 

lo S m dcr (X) dcr H> 

E^,{z,w;9) 

and consider the problem of minimizing f$[a](0) over all cr. Define 

(3.7) K$(0)= inf &Jcr](0); 

this quantity is related to M$(0) via (|3.5I) . in view of Theorem 1.3 p. 145]. 
Any a G P C (C) which achieves the maximum in (|3.7|) is called an extremal (or 
equilibrium) measure. In view of 13. 4f) . which follows from (|1.1|) . it follows from 
Theorem 1.3 of |191 p. 27] that there exists a unique extremal measure, which 
we will denote by a, or, when it necessary to indicate the dependence on the 
parameter 9 as well as on the weight $, by a® t g. 



4. The extremal measure and Hele-Shaw flow 

A change of parameters. For reasons of convenience, we will now use r = 1/(26*) 
in place of 9 as a positive real parameter. As before, $ : C — > K is C 2 -smooth, and 
subject to the growth requirement Hl.ljl for every choice of the real parameter A. In 
particular, l|3.4(l holds with 9 = l/(2r). As we mentioned in the previous section, 
by Theorem 1.3 ^3 p. 27], there exists a unique extremal measure a € V C (C) which 
achieves the infimum in (|3.7|) . For <r, the associated logarithmic potential 




zeC, 
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is locally bounded in C (Theorem 4.3 ^] p. 51]). So, for instance, the extremal 
measure a has no point masses. To get an understanding of this extremal measure, 
we make a detour to obstacle problems and Hele-Shaw flow. A general reference 
for Hele-Shaw flow is the recent book of Gustafsson and Vasiliev [JJ. 

An obstacle problem. Consider the following function: 

V(z) = -$(■?), zeC. 

For < t < +oo, let SH T (C) denote the collection of (extended real- valued) 
functions that are superharmonic in the whole plane C, and which are of the form 

-2t log I z I +R{z), 

where R(z) is superharmonic in a punctured neighborhood of infinity, and bounded 
from below there. If we have two functions /i, fi G SH T (C), then the minimum of 
the two - min{/i,/2} - is in SH r (C) as well. In addition, the decreasing limit of 
functions in SH r (C) remains in SH T (C) unless it degenerates to — oo identically. 
We define V T to be the least majorant (or lower envelope) ofV in the class SH T (C). 
In other words, V is an obstacle, and V T solves the obstacle problem. 

The extremal measure and the obstacle problem. The next proposition explains 
the relationship between the obstacle problem and the solution da to the opti- 
mization problem l|3.7[) with 9 = l/(2r). For convenience of notation, we shall 
write <7$ iT in place of o : * ) i/(2 T ) when we need to indicate the dependence on the 
parameter r. 

Proposition 4.1. The logarithmic potential for the extremal measure a = a$ !r is 

L[*](z) = ±-V T (z) + C T , 

where C T is the constant 

C - = ^r\ f // log—^da(z)dcr(w) 

2t{Jc JJc 2 \z-w\ 



Proof. This follows from Theorem 4.1 [03 p. 49]. □ 

We recall that the function V = — $ is assumed to be of class C 2 throughout 
C. We shall need the class C 1,1 , which consists of continuously differentiable func- 
tions whose gradient is locally Lipschitz. In other words, all second order partial 
derivatives of the function are locally bounded. 

Proposition 4.2. The envelope V T is of class C 1 ' throughout C. 

Proof. In view of Proposition l4~Tl the envelope function V T is harmonic off supp a. 
Since the extremal measure a has compact support, it follows that we can find 
an open circular disk D(0, R) of radius R, < R < +oo, around the origin, such 
that supp a is compactly contained in ED(0, R). Along the boundary T(0, R) of the 
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disk, V T is harmonic, and in particular smooth. We consider the following obstacle 
problem. Suppose v is superharmonic on O(0, R) and C^-smooth on 0(0, i?), and 

f v(z) > V(z), z e B(0,i?), 



(4.1) 



v(z) = V T (z), zeT(0,R). 



The lower envelope of all such v is denoted by V T . A simple argument shows that 
V T coincides with V T . Indeed, we easily see that V T < V T on B(0, R), with equality 
on T(0, R); moreover, if we extend V T to all of C by declaring that it should equal 
V T off B(0, R), we get a superharmonic function, which belongs to SH T (C). But 
- by definition - V T is the smallest such majorant, and we must have equality: 

V T = V T . 

In view of the known smoothness properties of solutions of obstacle problems 
of the type 1)4. lfl . it follows from the C 2 -smoothness of V that V T is of class C 1 ' 
on B(0, R). In the rest of the plane, it is harmonic, and therefore automatically of 
class C 1 ' 1 . The basic references on obstacle problems are Chapter 1 in [S], or the 
paper by Caffarelli and Kinder lehrer. 

The proof is complete. □ 

We introduce the notation 
(4.2) V< s >(t) = {z€C:V(z)<V t (z)X; 

this is an open set which contains a punctured neighborhood of the point at infinity. 

Proposition 4.3. We have 

{zeC: A$(z) < 0} C £>#(r). 

Proof. Let zq S C be a point with AQ(zq) < 0, and suppose that zq ^ 2?$(r). 
Then 

V(z ) = V(z Q ) = -$(z ), 

so that if we put U = V — V, we get U (zo) = 0, while U > everywhere, and U is 
strictly superharmonic in a neighborhood of zq. By the strong maximum principle, 
then, this is not possible. The result is immediate. □ 

We also need the "harmonicity" set of 

H^{zeC: A$(z) = 0}, 

which is a closed set in C. If X is a Borel measurable subset of the plane C, let 
us agree to say that z 6 C is a pseudo-interior point for X if there exists a small 
open disk D around zq such that 

\Dnx\ 2 = \D\ 2 

holds, where | ■ 1 2 is the operation of taking the Lebesgue are measure of the given 
set. For instance, all interior points of the set X are pseudo-interior points of X. 

We are now in a position to reformulate in precise terms our second theorem 
from the introduction fTheorem ll.2|) . 
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Theorem 4.4. In the sense of measures, we have 

dd$, T {z) = ^— lc\u*(r)(^) A$(z) d\ 2 (z), z E C. 

In particular, the support of the measure a$ iT consists of all points ofC which are 
not pseudo-interior for the set 2?$(r) U 7i$. 

Proof. We get from Proposition ^. II that 

(4.3) -27rdff(z) = Ai[a](z)dA 2 (z) = ^ AVJz) dX 2 (z). 

2r 

According to p. 53], we have 

(4.4) AV T {z) = AV(z) = -A$(z), z £ C\T>${t), 

in the almost-everywhere sense. On the non-coincidence set, however, V T is har- 
monic: 

(4.5) AV T (z)=0, zeDt(r). 

This follows by a standard Perron process argument (see, for instance, Proposition 
2.2 in 22]). By Proposition ^. 21 the function AV T is in L^ C (C), so that we see from 
(PI) . E3I), and gSJ) that 

da<s, tT (z) = - — 1c\d*(t)(2) A$(z), z e C, 

as asserted. The statement regarding the support of <r$ )T is an easy consequence 
of this identity. The proof is complete. □ 

The obstacle problem considered here share many features in common with 
the Hcle-Shaw flow domains considered by Hedenmalm and Shimorin in |12j How- 
ever, the time parameter t flows in the opposite direction as compared with the 
time parameter t used in |12j . This means that any cusps of the compact set 
C \ 2?$(r) should point outward. Also, if $ is real-analytic, the smoothness analy- 
sis of based on Sakai's work [2U| done in ^5] carries over to the situation treated 
here. 

The domains T>(t) grow as r decreases, and the "harmonic moments" are 
preserved, with the exception of the first; this is the content of the following 
proposition. 

Proposition 4.5. Let < r < t' < +oo. Suppose h is harmonic and bounded in 
2?$(t) ; with an extension to C that is locally of Sobolev class W 2 ' 2 . We then have 
the equality 

4tt (t' - r) h(oo) = [ h(z) A<S>{z) dX 2 {z). 



Proof. This follows from a suitable application of Green's formula, analogous to 
what is done in D 
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If the domains D$(r) are smooth and vary smoothly with r, then, by argu- 
ing as in PP- 188-189, we see that the assertion of Proposition 14.51 has the 
interpretation that the boundary dT>§ (r) propagates with velocity proportional to 
a weight times the gradient of the Green function with singularity at infinity. 

The growth of the complementary sets C \ T>$, (r) as r increases is quite an 
interesting process. For r close to 0, we should expect each such complementary 
set to be localized close to the points where the minimum of the potential function 
$ is attained. If there is only one minimum point zo, and the function $ is convex 
in a neighborhood of zq, then we can show that a small neighborhood of zo is 
contained in C \ D$(r), for each fixed r, by comparing with a concave majorant. 
An interesting question seems to be whether 

n z>*(t)=0. 

0<T< + OO 

It is easy to see that a necessary condition for this to happen is that <& be subhar- 
monic everywhere. We do not know to what extent this condition is sufficient. 



5. Condensation of quantum Hele-Shaw flow 

The continuous limit of eigenvalues of normal matrices. For positive real (3, let 
/i/3,4> denote the measure 

(5.1) dfip^(z) = exp { - /3®{z)} d\ 2 (z), zeC, 

where we recall that dA2 is area measure in the plane. Throughout this section, 
we shall use - as before - the following scaling choice of (3 as we increase N: 

(5.2) p = 2 e(N-l) = ^^, 

T 

where r is a fixed positive real parameter. Also, we recall the previously used 
notation 

Z[l, n ] = ( z l, ■ ■ ■ , z n), Z [1,N] = ( z l) ■ ■ ■ , z n)- 

in) 

Given this setup, we consider the n-point correlation measure Y ^ N with /i = /i/3,$ 
and P given by l|5.2|) . We want to understand the asymptotic behavior of this 
measure as JV — > +oo. This means that we should understand the impact of 
working with the L 2 norm in place of the L°° norm in the setting of l|3.1(l . 

We recall that $ is a C 2 -smooth real- valued function, which tends to infinity 
at infinity at a pace prescribed by (|l.lfl . It follows from this that there exists a 
compact subset K of C such that 

(5.3) K$(e) + 1< log ' 



E$(z, w; 6) 



9 \${z) + $(id)1 + log , r , (z, w)eC 2 \(Kx K), 

\z — w\ 
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where we recall the definition 13.7|l of k$ (9) . In analogy with the energy function 
(13.61) . we introduce the function S\, [zri.jvi] ($)> 



(5-4) 41^,]^) = ^^ J[ - 



{zj,zk;0) 



where E^(z, w; 9) is given by (|3.3|) . and 9 = l/(2r). We recall that 

(5-5) l ° g AU N (9) ~ £ *hwW)> Z H^] e C " 

with equality only at the weighted Fekete points. For positive real e, let 

2kiv(e, 0) - jz^jv] G : 4[a [llJV] ](0) < log + 4- 

We need to know that the proportion of points in 2l$jv(e, 0) which stays in the 
compact set K converges to 1 as N — > +oo and e — ► 0. 

Proposition 5.1. Suppose Zhjv] G 2ljv(e, #)> < e < 4. Le£ TVr- denote 

the number of indices j for which Zj G if. Then, for sufficiently large N , we have 

**>l-2e. 

AT ~ 



Proof. Let X C {1, . . . , TV} be the subset of all indices j for which Zj € if, and 
let y be the complement in {1, . . . , N}. We split the sum defining £|, [%[i,n]\ (&) 
accordingly: 



2 i 1 



1 V- j 1 



iV(iV 

The sum of the last two terms is estimated from below by 

1 N K (N K -1) \ ( , 
1 N(N-l) JM^ + 1 )' 
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while the first term of the right hand side is estimated in the following manner: 
1 



log 



AL 



<S>,N K 



We now see that for Z[i,at] € 2ljv(e,#), 



N(N-l) M^ Nk (9) V N(N-1) 



< £l{z [hN] }(6) <log- ' 



It follows that 

(5.6) l + K g(fl)-log — -i— -e < ~. 1) (l + „*(0)-log — — , 

M<s> tN {0) N[N-1) \ M$ tNK (9)J 

For large values of AT, the left hand side is positive, and hence the right hand side 
is, too. We quickly rule out the possibility that the first on the right hand side 
is negative (which corresponds to Nk = 0), in which case we are left with both 
factors on the right hand side being positive. We get from l|5.6|) that 

(5-7) iV ±"" 1) >l- 



N(N — 1) " l + K # (0)-log 



if we recall that M$>^(9) < M^, i n k (9). Next, we note that Nk — > +oo as N tends 
to infinity; therefore, for sufficiently large N, 

The assertion of the proposition now follows from (|5.7|) . □ 

Lemma 5.2. Let a £ 'Pc(C) 6e absolutely continuous with respect to two-dimensional 
Lebesgue measure, 

da{z) = S(z)d\ 2 (z), 

where S > is area-summable and S log + 5 is area-summable as well. The nor- 
malization constant Zjy in the definition o/dIIjv,^ J $(z[ li jv]) then has the following 
bound: 

Z N >expi-N(N-l)£<s>[o](9)-N J log S(z) da(z)j, 
with the understanding that logS'dcr vanishes off the support of a. 
Proof. We recall first the definition of Zpj: 

(5.8) Z N = ^CXpj -/3^$(^)| |A(z [1 , A r ] )| 2 dA 2A r(z [liA r ] ). 
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Let S denote the set where S(z) > 0; as a has compact support, we may assume 
that the set £ is bounded. We then have 



Z N > ( exp( -ffVsOj)) |A( Z[1;JV] )| 2 1 — -dcr(zi)---dcxOjv) 

>ex P { / f-/?V$(z,) + 21og|A(z [1 . Ar] )| - Vlog5(z 

x da(zi) ■ ■ ■ dcr(zAr)|' , 

where the second estimate is due to Jensen's inequality. In view of (I5.2|l and the 
definition of the energy functional, we have 

Z N > exp j - N(N - 1) £ # [<x] (6>) — iV / log dcr(z) 1 , 



as claimed. □ 

Remark 5.3. We note that in view of Theorem 14.41 the extremal measure a = a 
enjoys the regularity property of Lemma 15.21 

The next proposition is central to the argument; its proofs mimics Johans- 
son's method for eigenvalues of self-adjoint matrices |13) . 

Proposition 5.4. For sufficiently large N, we have 



Proof. For z [1)Ar] £ C \ a w (e, 0), we have 

> k$(0) + e. 

We recall the estimate (ll.lf) from below of the weight <I>, which for an appropriate 
choice of the parameter A gives 



1 N 
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for a positive constant C that only depends on 9 and <I>. We form a convex com- 
bination of these two estimates (0 < 7 < 1): 

(5.9) 4 [z lhN] ] (9) > (1 - 7) [«* (0) + e] - C 7 

jv 

+ ^^log(i + |z,| 2 ), ZM eC\%M). 
j=i 

We recall the formula defining IIjv,/3,#, while keeping in mind the identity 13. 2[) . 

1 



(5.10) dIIjv, J 3,$(z[i.iv]) = ~ 



dA 2 Ar(z[ 1 ^ A r]) 



= J- cxp{-7V(7V-l)4[z [l!jv] ](0)}dA 27V (z [liJV] ). 
So, on the set \ 2ljy(e, 9), we get, in view of (|5.9|) . 

dn iv ^,$(z [1 , iv] ) < ^ cxp j -JV(JV-l)(l-7)[/s«(d)+e] 



JV 1 

+ C7V(iV - 1)7 - 7 (JV - 1) J] log (1 + | Zj | 2 ) A\ 2N (z [1>N]j 

3=1 J 

By Lemma 15. 21 we have 

Z N > expj - N(N- 1)k$(0) -Bn\, 



where B is the real number 

f 9 
B= logA$(z)dCT(z) + log-. 

JC T 

As we combine this with (|5.9|l . we arrive at 

dU N , ^(z [ltN] ) <e(xp|-JV(JV-l)[7«»(fl) + (l-7)e] +BN 

N v 

+ C 7 jV(iV - 1) - 7 (JV - 1) £ log (1 + \ Z j | 2 ) [ d\ 2N (z [hN] ) 
= exp I - jV(jV - 1) [7 k$(0) + (1 - 7 )e] + £ TV + C7jV(jV - 1) 

JV 

on the set C N \ 2ljv(e, Considering that 

(l + N 2 )-^dA 2 (z)- 11 



7(JV- 1) - 1' 
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provided that 1/(N — 1) < 7 < 1, we find that 

n Njl3 ^(C N \ % N (e, 6)) < exp ( - N(N - 1) [7 k*(0) + (1 - j)e - C7] + BN 



n 



^(N- 1) - 1 

The dominant contribution in the expression which is exponentiated is 

-N(N - 1) [7 k^(9) + (1 - 7 )e - C 7 ] . 

We would like to pick 7, 1/(A — 1) < 7 < 1, such that 

0<7«$(6') + (l-7)e-C , 7. 

This is possible; indeed, without loss of generality, we may assume that C is greater 
than k$(0), in which case one such choice is 

e 

7 = 2(C-k*(0)+c)' 
This value of 7 yields the assertion of the proposition. The proof is complete. □ 

For a point zrj^ £ C , we define the associated weighted sum of point 
masses ct[z[ 1JV ]] £ P C (C) by the formula 



(5.11) da[z [1>iV] ](z) = -i^d^( 2: ) ! zeC, 

i=i 

where 5 tt means the Dirac point mass at w G C. Also, let Cf,(C) = C(C) n L°°(C) 
denote the space bounded complex-valued continuous functions on C. 

Proposition 5.5. Suppose o~n — o~\z[i N ^ is as above, with 

Z [1,JV] = ( z i, ■■■,zn) e < C N - 

Suppose, moreover, that 

4[*[i,JV]](0)-/«*(0) 
as A — > +00. TTien <7/v converges to a, the extremal measure, in the weak-star 
topology, as N — > +00. /n o£/ier words, for each f s Cf,(C), we Ziaue 



/(z) d<7 n{z) — 51 / /(^) df?(z) as A — > +00. 



Proof. First, fix a small but positive e. By assumption, for sufficiently big A, 

Z[i,N] G Stjv(e,0), 

so that by Proposition l5.ll <jn(K) > 1 — 2e. As e can be made as small as we like, 
and each o~n is a probability measure, we see that 

ctn(K) -> 1 as A -^+00. 
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The space of all finite complex Borel measures on K is a Banach space, with 
weak-star compact unit ball. This means that each subsequence of the sequence 
0i|k> g^Ki Gs\K) ■ ■ ■ has a weak-star convergent subsequence. We shall show that 
any such limit coincides with a, from which the assertion follows. 

Without loss of generality, then, by passing to a subsequence, we may assume 
that the sequence <jn\k, N = 1,2,3, converges weak-star itself to a limit, which 
we call cr. The weak-star limit of <jn is then also a, given that ctjv(C \ K) — > 
as N — > +oo. By testing with the function / = 1, we see that cr is a probability 
measure supported on K, so that a €E V C (C). We claim that £$[<?] (0) < k&(6). 
Once this is established, the equality a = a follows from the uniqueness of the 
extremal measure. We use a cut-off argument: we note that 

2 Q N 1 i 

so that if L is a real parameter, 

90 N i N r 1 1 



N -1 



Now, let N — > +oo, so that we get, in view of our assumptions, 

k< s ,(6)>29 / $(z)dcT(z) + / mini log- — - — r , L \ da(z) da(w). 
Jc Jc* I \z-w\ J 

As we let L tend to +oo, the right hand side approaches £$[<?] (#), whence the 
claim is immediate. □ 

As in the introduction, we let n^jg 4 be the n-point marginal distribution 

(n) 

measure and T N '» $ the n-point correlation measure for the probability distribu- 
tion measure Un,/3,$- We keep f3 connected with r, N, and via l|5.2[) . and fix 
(or, if you like, r) while A grows. 

Let Cb(C") = C(C")nL 00 (C rl ) denote the Banach space of bounded complex- 
valued continuous functions on C™ . We arrive at a precise reformulation of our first 
main result (Theorem II. 1(1 . 



Theorem 5.6. We have 
while, as N —> +oo, 



<k*(C") = 1, 



Moreover, for each f £ Cb(C"), we have, as A 

T («) 



( 5 - 12 ) / /( Z [l,«]) dn 7V,/3,<f ( Z [l,"]) ^ / / ( z [l,n] ) da(zi ) • • • dc ? (z„) 
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In other words, in the weak-star topology of measures, we have, as N — > +00, 

dn 7V,i3,<£>( z [i,«]) dCT(zi) ■ ■■da(z„). 
Here, a — a<& t g is the extremal measure. 

Proof. We note that the total mass of the measure nj^L # is 1. 

As a second step, we prove (I5.12f) under the slightly more restrictive assump- 
tion / G C C (C), which means that the test function / has compact support. 

Let ? be a permutation of {1, 2,3,..., N}. Then, due to the symmetry prop- 
erties of Hn,0,$, 

/ f( z [l,n]) dn 7V/3 *( Z [1,«]) = / ■ ■ ■ dnAr >/3i $(z[ liAr ]) 

. . . , Z c („)) dIl7v,/3,*(z[l,Ar]), 



from which we quickly deduce that 

/(z[l,n]) dn^ |$ (z[l,„]) = J CN ■ ■ ' ' z ?(«)) dn JV,/3,*(z[l,iV]), 

where the sum runs over all permutations of {1, 2, 3, ... , N}. We fix a small positive 
e, and split the integral: 

(5.13) / f( Z[1 , n] )dn ( ^( Z[1 , n] ) 

Jc 

+ aTiE / /( z 5 (i)v,z 5 w)<in w (z[i,jv]). 

By Proposition [Q] the last term is bounded in modulus by 

(5-14) e-^-^ll/IU^Q, 

for large N. To deal with the first term on the right hand side, we should under- 
stand the behavior of 

(5-15) ^y^/(z ?(1) ,...,z ?( „)), i [1)N] €% N (e,6). 



Let us consider the simplest case n = 1 first. Then (|5.15() amounts to 

JV 

A' 



1 N 

(5-16) ^E/fe)' Me, 0). 



j'=i 
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By letting e approach slowly as TV — > +00, we may ensure that l|5.14|l tends to 
as TV — > +00, while l|5.16|l approaches 

Jc 

It is easy to check that the latter statement entails 

1 f N f 

Jj\ / y]/( z 3) dII -iV^ 1 *(z[i 1 jv]) -> / /(z)da(z) 

iV! Jsi N (e,6) ~{ JC 

as TV — > +00, if we allow e to approach slowly. 

The remaining case n > 1 is handled in an analogous manner: Proposition 
15.51 should be replaced by a multidimensional analogue, which we get by iterated 
integration. It is useful to keep in mind that for large TV and fixed n, the following 
collections of sequences are asymptotically (as TV — > +00) the same: 

(f(l),...,c(n)) and . . . ,j(n)), 

where ? runs over all permutations of {1, . . . , TV}, and j runs over all functions 
j : {1, . . . , n} — > {1, . . . , TV}. For instance, 

N n (N-n)\ 

— ; ► 1 as iV -> +00. 

iV! 

There are a couple of assertion that remain to check. If we notice that the 
support of a is actually contained in the interior of the compact set K, then, 
by choosing a smooth cut-off function / = x S C C (C") with < x{ z [i,n]) < 1 
everywhere, x( z [i,n]) = 1 011 ( su PP^)"i an d xC- 2 ) = off if™, we find from l|5.12() 
that 

hminf flW (K») > 1. 
Since we are dealing with probability measures, it follows that 

hm U ( ^(K n ) = l- 

as claimed. 

It remains to verify that 15.12|) holds for all / S Cb(C"). To this end, we 
use a smooth cut-off function x similar to what we defined above, and write / = 
X / + (1 — x) /■ Then %f € C C (C), while the integral of the remainder (1 — x) / 
tends to as TV — > +00, by a simple estimate. 

The proof is complete. □ 

Remark 5.7. Let us think of the TV-tuple znjn = {z\, ... , z/v) as a random variable 
taking values in C N , with distribution measure IIjv,/3,$. The marginal distribution 
of the first n (1 < n < TV) coordinates Z[x >n ] = ( z i? ■ ■ ■ 1 z n) is then given by the 

(n) 

measure IT^ ^ ^; to stress the dependence on TV, let us write 



5 [l,n] I AT — (Zi\n, ■ ■ ■ : Zti\n) 
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for the marginal random variable. Now, if (3 grows with N according to l|5.2[) . 
for some fixed positive 6, Theorem 15.61 savs that as N — > +oo, Zi\n, ■ ■ ■ , Zn\N 
become identically distributed independent random variables, each having a$ t g as 
distribution measure. In particular, the covariance type matrix (with fi = /i/3,$) 

K M,Ar( Z [Ml]: Z [l,«]) = ( K H,N(Zj,Zk)) 

becomes asymptotically diagonal as N — > +oo. 

Example. We need to give an example to illustrate the result. We consider the 
rather trivial case of &(z) = \z\ 2 , and put n = 1. Then A<j>(z) = 4 is constant. 
A computation shows that with fi = /i/3.$ given by (|1.2() . the reproducing kernel 
function for the polynomial subspace is 

K MtN (z,w) = - 2^ q — • 



1 — ' 7 



We then have 



dn&W« =^^fe.)e-^l 2 dA 2 (z) = ^e-^l 2 £ M^dA 2 (z), 

77 i=o J ' 

while it is easy to check that 

X>*(t) = G C : |z| > 

In view of this, the content of Theorem 15.61 in this simple case is that, for fixed 
positive t, the function 

. (J v-D N Vr V I / ^-1)M 2NJ ' 



which is real- valued, with values between and 1, tends to 1 as ./V — > +00 for 
\z\ < y/r, and to for \z\ > y/r. This fact is of course well-known. In 0, the 
domain 2?$(t) is computed explicitly and shown to be an ellipse in the more 
general case (still with constant A$(z) = 4) 

$(z) = \z\ 2 + aRe{z 2 ), -1 < a < 1. 
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